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ABSTRACT 

We announce the discovery of 1 - 100 GeV gamma-ray emission from the archetypal TeV pulsar wind 
nebula HESS J1825— 137 using 20 months of survey data from the Fermi Large Area Telescope (LAT). The 
gamma-ray emission detected by the LAT is significantly spatially extended, with a best-fit rms extension of 
g = 0.56° ± 0.07° for an assumed Gaussian model. The 1 GeV - 100 GeV LAT spectrum of this source is 
well described by a power-law with a spectral index of 1.38 ± 0.12 ± 0.16 and an integral flux above 1 GeV 
of (6.50 ± 0.21 ± 3.90) x 10~ 9 cm -2 s — 1 , The first errors represent the statistical errors on the fit p arameters, 
while the second ones are the systematic uncertainties. Detailed morphological and spectral analyses bring new 
constraints on the energetics and magnetic field of the pulsar wind nebula system. The spatial extent and hard 
spectrum of the GeV emission are consistent with the picture of an inverse Compton origin of the GeV- TeV 
emission in a cooling-limited nebula powered by the pulsar PSR J 1826— 1334. 

Subject headings: pulsar, pulsar wind nebula, ISM: individual objects (G18.0-0.7, HESS J1825— 137), pulsars: 
individual (PSR B 1823 -13, PSR J 1826 -1334) 



1. INTRODUCTION 

Pulsars dissipate their energy via magnetized particle winds 
(consisting of electron/positron pairs) with the confinement of 
this particle wind-outflow leading to the phenomenon of pul- 
sar wind nebulae (PWNe). PWNe develop when the parti- 
cle wind collides with its surroundings, in particular with the 
slowly expanding supernova ejecta, and form a termination 
shock. Even though many PWNe are bright enough to be re- 
solvable both spatially and spectrally at wavelengths from ra- 
dio to very high energy (VHE; E > 100 GeV) gamma-rays, 
many open questions remain. How is the rotational energy of 
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the pulsar converted into a relativistic particle wind? What is 
the bulk Lorentz factor of the wind, what is the mechanism by 
which the particles are accelerated at the termination shock? 
What is the partition between magnetic and particle energy in 
the wind? To address these questions, PWNe can be seen as 
laboratories for relativistic astrophysics with the advantage of 
having a well-localized (and usually well-characterized) cen- 
tral energy-source, the spin-down power of a neutron star. 

Many studies of PWNe have been conducted in the X- 
ray b and with instruments suc h as XMM-Newton or Chan- 
dra dKargaltsev & Pavlovll2008l) . but recent gamma-ray ob- 
servations with H.E.S.S., MAGIC and VERITAS have con- 
tributed significantly to the understanding of PWNe. To- 
gether with the assumption that the X-ray emission is syn- 
chrotron whereas the higher energies are Inverse Compton 
(IC) emission, these observations allow us to derive the 
magnetic field strength in the nebula. For the PWNe seen 
in TeV gamma-rays, typical magnetic field strengths are ~ 
10 /iG (Ide Jager & Djannati-Ataf 20091) . In turn, the lifetime 
of the synchrotron-emitting electrons can be determined from 
the magnetic field strengths. ForkeV X-ray emitting electrons 
the lifetime of less than 100 years is significantly shorter than 
the age of the central pulsar ( 10 3 — 10 5 years) whereas the TeV 
emitting electrons have lifetimes of order 10 4 years, compara- 
ble to the pulsar age. This has the consequence that the X-ray 
PWNe should be spatially very compact, whereas the PWNe 
as seen in TeV should be significantly larger. The discov- 
ery of energy-dependent morphology at TeV energies in the 
archetypal system HESS J1825— 137 (surrounding the ener- 
getic pulsar PSR J 1 826— 1 334) confirmed this basic picture by 
demonstrating that the emission is dominated by 'relic' elec- 
trons from the earlier epochs of the nebula in which the pulsar 
was spinnin g down more rapidly, r eleasing more energy into 
the system dAharonian et al.ll2006l) . The gamma-ray studies 
of PWNe also showed that the overall efficiency of conver- 
sion of rotational energy into TeV emission in these systems 
is extremely high (^> 10% in some cases). 

The recently-launched Large Area Telescope (LAT) on 
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board the Fermi Gamma-ray Space Telescope has established 
pulsars as the most numerous c lass of identified Galactic 
GeV gamma-ray emitting objects ( Abdo et al. 2010b) and has 
also shown that up to 10% of the rotational power from pul- 
sars is emitted in the Fermi-LAT energy band. The spec- 
trum of the pulsed comp onent cuts off in the sub-10-GeV 
range (Abdo et al. 2010a). As detailed above, the spectra of 
PWNe extend out to TeV energies making PWNe the most 
numerous Galactic objects at those energies. While for a 
few TeV sources there is a clear morphological match (e.g. 
with the X-ray emission) that leads to the identification with 
a PWN, for most TeV sources the situation is less obvious. 
Generally, a positional match with a pulsar (along with the 
requirement that the energetics of the pulsar is sufficient to 
power the TeV emission) is used to suggest an association. 
The situation is, however, further complicated by the fact that 
for many of the TeV-detected PWNe, in particular for the ones 
that are associated with middle-aged rather than young pul- 
sars, the gamma-ray emission is offset from the pulsar. 

HESS J 1825 — 137 is one of the prime examples for such 
an offset PWN system. The system seems to be powered by 
the energetic radio pulsar PSR J1826 — 1334. High resolution 
X-ray observations of the PWN with XMM-Newton showed 
a compact core with a hard photon index (T = 2) °f 

size 30" embedded in a larger diffuse structure of extension 
~ 5' extending to th e south of the pulsar with a softer photon 
index of V ~ 2.3 dGaensler et al.ll2003l) . The TeV gamma- 
ray emission has a much larger extent (~ 0.5°) and shows 
a similar softening of the photon index from T = 2.0 close 
to the puls ar to T = 2.5 at a dis tance of 1° to the south of 
the pulsar ( Ahar onian et al.ll2006l) . Based on hydrodynamical 
simulations (Blondin et al. 2001) the asymmetric nature of the 
emission can be explained by dense interstellar material to 
the north of the pulsar. The reverse shock of the supernova 
explosion from that direction might have been pushed inward, 
interacting therefore relatively early with the PWN and thus 
pushing the X-ray and TeV gamma-ray emission mainly to 
the south. 

Here we report on GeV gamma-ray observations of the 
HESS J1825-137/PSR J 1826- 1334 system with Fermi- 
LAT. 

2. LAT DESCRIPTION AND OBSERVATIONS 

The LAT is a gamma-ray telescope that detects photons 
by conversion into electron-positron pairs and operates in 
the energy range between 20 MeV and more than 300 GeV. 
It is made of a high-resolution converter tracker (for direc- 
tion measurement of the incident gamma-rays), a CsI(Tl) 
crystal calorimeter (for energy measurement) and an anti- 
coincide nce detector to ident ify the background of charged 
particles (lAtwoo d et al. 2009). In comparison to EGRET, the 
LAT has a larger effective area (~ 8000 cm 2 on-axis above 
1 GeV), a broader field of view (~ 2.4 sr) and a superior an- 
gular resolution (~ 0.6° 68% containment at 1 GeV for events 
converting in the front section of the tracker ). Details of 
the ins trument and data processing are given in Atwo od et al. l 
(2009). The on-orbit calibration is described in Abd o et al.l 

The following analysis was performed using 20 months of 
data collected starting August 4, 2008, and extending until 
April 21, 2010. Only gamma-rays in the Pass 6 Diffuse class 
events were selected (i.e. with the tightest background re- 
jection), and from this sample, we excluded those coming 
from a zenith angle larger than 105° to the detector axis be- 



cause of the possible contamination from secondary g amma- 
rays from the Earth's atmosphere (lAbdo et al.l l20Q9bb . We 
have used P6_V3 post-launch instrument response functions 
(IRFs), which take into account pile-up and accidental coin- 
cidence effects in the detector subsystems. 

3. TIMING ANALYSIS OF THE PULSAR PSR J1826— 1334 

The pulsar PSR J1826-1334 (also known as 
PSR B1823— 13) was discovered in the survey of 
IClifton et aT] d!992l) . Its spin period of 101.48 ms, character- 
istic age of 21 kyrand spin-down power of 2.8 x 10 36 ergs _1 
are very similar to those of the Vela pulsar. The distance 
to the pulsar a s derived from the d ispersion measure is 
(3.9 ± 0.4) kpc dCordes & Lazioll2002l) . Although the pulsar 
was not reported in the Fermi-LA T catalog of gamm a-ray 
pulsars with 6 months of data (Ab do et all l2010al) . we 
performed a temporal analysis on this increased Fermi-LAT 
dataset using a new timing solution. A total of 162 obser- 
vations of P SR J1826-1334 were m ade at 1. 4 GHz using 
the P arkes dWeltevrede et al.1 120091). Lo vell dHobbs et all 
|2004|) and Nancav dCognard et al.l 1201 Ol) radio t elescopes. 
The TEMP02 timing package dHobbs et al J 120061) was then 
used to build the timing solution. We fit the radio times 
of arrival (TOAs) to the pulsar rotation frequency and first 
three derivatives. We whitened the timing noise with four 
harmonically related sinusoids, using the "FITWAVES" 
functionality of the TEMP02 package. The post-fit rms is 
292.4 /is, or 0.2 % of the pulsar phase. This timing solution 
will be made available through the Fermi Science Support 
Center 14 . 

Photons with an angle 8 < max(5.12° x 
(£yi00MeV)~°' 8 ,0.2°), where E is the energy of the pho- 
ton, from the radio pulsar position, R.A. = 276.55490°, decl. 
= -13.57967° (J2000), were selected and phase-folded us- 
ing the above mentioned ephemeris. The energy-dependence 
of the integration radius is a satisfactory approximation of the 
shape of the LAT Point Spread Function (PSF), especially at 
low energies. 

The H-Test values, as defined in de Jageret al. ( 1989) and 
obtained from the analysis of the pulsed emission, correspond 
to a significance well below 2er for each tested energy band 
(30 MeV - 300 GeV, 30 MeV - 100 MeV, 100 MeV - 300 
MeV, 300 MeV - 1 GeV, > 1 GeV). No significant pulsation 
is detected with the current statistics. We fitted a point source 
at the position of the pulsar PSR J 1826 — 1334 and derived an 
upper limit between 100 MeV and 1 GeV of — 3.1 x 10~ 8 
cm~ 2 s -1 , well below typical gam ma-ray fluxes repo rted for 
pulsars detected by Fermi-LAT dAbdo et alj l2010al) . This 
ensures that any emission from the PWN HESS J 1825- 137 
will not be contaminated by pulsed gamma-ray photons from 
PSR J1826- 1334. 

4. ANALYSIS OF THE PULSAR WIND NEBULA HESS J1825- 137 

The spatial and spectral analysis of the gamma-ray emis- 
sion was performed using two different methods, gtlike 
and Sourcelike, gt like is the maximum-likelihood method 
(Matt oxet alj [l996) implemented in the Fermi SSC science 
tools. This tool fits a source model to the data along with 
models for the instrumental, extragalactic and Galactic back- 
grounds. In the following spectral analysis, the Galactic 
diffuse emission is modeled using the ring-hybrid model 

14 FSSC: http://fermi.gsfc .nasa.gov/ssc/data/access/lat/ephems/ 
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gllJem.v02.fit. The instrumental background and the extra- 
galactic radiation are described by a single isotropic compo- 
nent with a spectral shape described by the tabulated model 
isotropic J.em_v02.txt. The models and their detailed descrip- 
tion are released by the LAT Collaboration I5 . Sources within 
10° of the pulsar PSR J1826- 1334 and found above the back- 
ground with a statistical significanc e larger than 5 a ar e ex- 
tracted from the source list given in lAbdo et ail (1201 Obi) , ex- 
cept for 1FGL J1825.7— 1410c which contributes partially to 
the gamma-ray emission of the PWN. A detailed analysis re- 
veals that the GeV emission is significantly extended com- 
pared with that of a point source. The morphological analy- 
sis was performed using Sourcelike, which is described in the 
Appendix. The extension test was done using a uniform disk 
and a Gaussian spatial model. Sourcelike can also be used to 
assess the Test Statistic (TS) value and to compute the spectra 
of both extended and point-like sources. The TS is defined as 
twice the difference between the log-likelihood L% obtained 
by fitting a source model plus the background model to the 
data, and the log-likelihood Lq obtained by fitting the back- 
ground model only, i.e TS = 2(£i — L ). 

4.1. Morphology 

To study the morphology of an extended source, a ma- 
jor requirement is to have the best possible angular resolu- 
tion. Therefore, we decided to restrict our LAT dataset to 
events with energies above 10 GeV. This also reduces the rel- 
ative contribution of the Galactic diffuse background. Fig- 
ure[T|(top) presents the LAT counts map of gamma-ray emis- 
sion around HESS J1825— 137, smoothed with a Gaussian of 
a = 0.35°. The image contains emission from the Galactic 
diffuse background and from nearby sources but in addition 
shows bright emission south of PSR J1826— 1334 coinciding 
generally with the region that is bright at TeV energies (de- 
noted by the H.E.S.S. flux contours). The Test Statistc (TS, as 
defined above) map above 10 GeV presented in Figure[T](bot- 
tom) supports this picture. This skymap contains the TS value 
for a point source at each map location, thus giving a measure 
of the statistical significance for the detection of a gamma- 
ray source in excess of the background. Clearly, significant 
emission to the south of PSR J 1 826— 1 334 is detected. 

We determined the source extension using Sourcelike with 
a uniform disk hypothesis and a Gaussian distribution (com- 
pared to the point-source hypothesis). The results of the ex- 
tension fits and the improvement of the TS when using spa- 
tially extended models are summarized in Table Q] The differ- 
ence in TS between the Gaussian distribution and the point- 
source hypothesis is TS C xt = 72 (which converts into a sig- 
nificance of ~ 8<7 for the source extension) for 10 GeV < E < 
100 GeV, which demonstrates that the source is significantly 
extended with respect to the LAT point spread function (PSF). 
The fit extension has a dispersion of a = 0.56° ± 0.07°. We 
support this conclusion in Figure [2] showing the background 
subtracted radial profile for the LAT data above 10 GeV (from 
the best source location determined for a Gaussian fit) and 
comparing this with the LAT PSF. Similar results are obtained 
assuming a uniform disk model. 

We have also examined the correlation of the gamma-ray 
emission with different source morphologies by using gtlike 
with assumed multi-frequency templates. For this exercise we 
compared the TS of the point source and Gaussian distribution 
parameters provided by Sourcelike with values derived when 

15 Fermi Science Support Center: http://fermi.gsfc.nasa.gov/ssc/ 



using the H.E.S.S. gamma-ray ex cess map as a morphological 
template ( Aharonian et al. 2006). The resulting TS values ob- 
tained from our maximum likelihood fitting are summarized 
in Table [2] Fitting a uniform disk to the data using the best 
location and size provided by Sourcelike instead of a point- 
source hypothesis results in TS ex t = 67, comparable to the 
improvement in TS between D and PS models in Table Q] 
Fitting a Gaussian model improves the TS by 80. Replacing 
the Gaussian with spatial template provided by the H.E.S.S. 
observations decreases the TS with respect to the Gaussian 
hypothesis (ATS = —42), implying that the LAT emission is 
not perfectly reproduced by the H.E.S.S. excess map (E > 
200 GeV). This is not completely surprising since H.E.S.S. 
reported an energy-dependent morphology for this source and 
Fermi is probing lower energy electrons. Thus, while the best 
match is with the Gaussian morphology, we cannot rule out a 
simple disk morphology. 

TABLE 1 

CENTROID AND EXTENSION FITS TO THE LAT DATA FOR 
HESS J 1825 — 137 USING Sourcelike FOR EVENTS WITH ENERGIES 
above 10 GeV. The difference in Test Statistic between a 
given spatial model and the point-source hypothesis is 

INDICATED BY TSext IN THE LAST COLUMN. 



Model 


Name 


in 


b(°) 


Radius (°) 


TS cx t 


Point Source 


PS 


17.62 


-0.82 






Disk 


D 


17.57 


-0.45 


0.67 ± 0.02 


72 


Gaussian 


G 


17.57 


-0.43 


0.56 ± 0.07 


72 



TABLE 2 

Comparison of model likelihood fitting results with gtlike 

FOR EVENTS WITH ENERGIES ABOVE 1 GEV. FOR EACH MODEL, WE 
GIVE THE NAME, AND THE TEST STATISTIC VALUE (TS). 



Model 


Name 


TS 


Point Source 


PS 


24 


Disk 


D 


91 


Gauss 


G 


104 


H.E.S.S. 




62 



4.2. Spectral analysis 

The following spectral analyses are performed using gtlike. 
The Fermi-LAT spectral points were obtained by dividing the 
1 - 100 GeV range into 6 logarithmically-spaced energy bins 
and performing a maximum likelihood spectral analysis in 
each interval, assuming a power-law shape for the source. 
For this analysis we used the Gaussian model from Table Q] 
to represent the gamma-ray emission observed by the LAT, 
as discussed in section 14.11 Assuming this spatial shape, the 
gamma-ray source observed by the LAT is detected with a TS 
of 104 (~ 10a) in the 1 - 100 GeV range. To determine the 
integrated gamma-ray flux we fit a power-law spectral model 
to the data in the energy range 1 GeV - 100 GeV with a max- 
imum likelihood analysis. This analysis is more reliable than 
a direct fit to the spectral points since it accounts for Pois- 
son statistics of the data. The spectrum of HESS J1825 — 137 
between 1 and 100 GeV, assuming the Gaussian model from 
Table Q] is presented in Figure [3] It is well described by a 
power-law with a spectral index of 1.38 ± 0.12 ± 0.16 and 
an integral flux above 1 GeV of (6.50 ± 0.21 ± 3.90)xl0" 9 
cm~ 2 s _1 . This is in agreement with results obtained inde- 
pendently using Sourcelike. The first error is statistical, while 
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FIG. 1. — Top: Fermi-LAT counts map above 10 GeV of the HESS J 1 825— 137 region with s ide-length 5°, binned in square pixels of side length 0.05°. 
The map is smoothed with a Gaussian of a = 0.35°. H.E.S.S. contours (Aharonian et al. 2006) are overlaid as gray solid lines. The position of the pulsar 
PSR J1826— 1334 and of the close-by 1FGL sources are indicated with green and white squares respectively. LS 5039 is visible in the South-East at position 
(RA, Dec) = (276.56°, —14.85°). Bottom: Fermi-LAT Test Statistic (TS) map for events with energy larger than 10 GeV on a region of 2.5° side length. The 
TS was evaluated by placing a point-source at the center of each pixel, Galactic diffuse emission and nearby sources being included in the background model. 



the second represents our estimate of systematic effects as 
discussed below and is dominated by the uncertainties on the 
Galactic diffuse emission in the 1-5 GeV energy range. With 
the current statistics, neither indication of a spectral cut-off at 
high energy nor significant emission below 1 GeV can be de- 
tected. 

Four different systematic uncertainties can affect the LAT 
flux estimation : uncertainties on the Galactic diffuse back- 
ground, on the morphology of the LAT source, on the effec- 
tive area and on the energy dispersion. The fourth one is rel- 
atively small (< 10%) and has been neglected in this study. 
The main systematic at low energy is due to the uncertainty 
in the Galactic diffuse emission since HESS J1825— 137 is 
located only 0.7° from the Galactic plane. Different ver- 
sions o f the Galactic diff use emission, generated by GAL- 
PROP (IStrong et all 12004). were used to estimate this error. 
The observed gamma-ray intensity of nearby source free re- 



gions on the galactic plane is compared with the intensity 
expected from the galactic diffuse models. The difference, 
namely the local d eparture from the b est fit diffuse model, is 
found to be < 6% (lAbdo et al.ll2010el) . By changing the nor- 
malization of the Galactic diffuse model artificially by ±6%, 
we estimate the systematic error on the integrated flux of 
the PWN to be 70% below 5 GeV, 34% between 5 and 10 
GeV, and <12% above 10 GeV. The second systematic is re- 
lated to the morphology of the LAT source. The fact that we 
do not know the true gamma-ray morphology introduces an- 
other source of error that becomes significant when the size 
of the source is larger than the PSF, i.e above 600 MeV for 
the case of HESS J1825— 137. Different spatial shapes have 
been used to estimate this systematic error: a disk, a Gaus- 
sian distribution and the H.E.S.S. template. Our estimate of 
this uncertainty is ^30% above 1 GeV. The third uncertainty, 
common to every source analyzed with the LAT data, is due 
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FIG. 2. — Background subtracted radial profile of the LAT data from the best-fit position provided by Sourcelike for a Gaussian source (1, b) = (17° .57, -0°.43) 
as reported in Table[T](E > 10 GeV). The best fit model, obtained for a Gaussian distribution, and the LAT PSF are overlaid as a blue solid and a red dashed lines 
respectively for comparison. The black dots represent the difference between the observed counts and the Galactic and extragalactic diffuse emission model. 
Nearby sources are not subtracted from the radial profile. 



to the uncertainties in the effective area. This systematic is 
estimated by using modified instrument response functions 
(IRFs) whose effective area bracket that of our nominal IRF. 
These 'biased' IRFs are defined by envelopes above and be- 
low the nominal dependence of the effective area with en- 
ergy by linearly connecting differences of (10%, 5%, 20%) 
at log(E) of (2, 2.75, 4) respectively. We combine these vari- 
ous errors in quadrature to obtain our best estimate of the total 
systematic error at each energy and propagate through to the 
fit model parameters. 

5. DISCUSSION 

A one-zone spectral energy distribution model provides a 
useful tool to investigate the global properties of the PWN. 
The nebula is non-uniform in the VHE regime (~ 1° at 1 
TeV and ~ 0.2° at 20 TeV) and p ossesses a bright ce ntral X- 
ray core observed with Chand ra dPavlov et al.ll2008|). XMM- 
Newto n dGaensler et al.l 12003) and Suzaku (lUchivama et al.l 
2009) which extends no more than 15' from the pulsar. The 
non-uniform X-ray and VHE morphologies likely stem from 
cooling losses by energetic electrons as they traverse the neb- 
ula, yet at lower energies in the uncooled regime the elec- 
tron spectral shape remains essentially constant with time, and 
hence also with position. As a result, even though a one-zone 
model cannot reproduce the energy-dependent morphology of 
the nebula, a one-zone model can nevertheless accurately re- 
produce the global flux from uncooled electrons. Electrons 
inverse-Compton scattering off the CMB require energies of 
~ 2 TeV in order to produce photons of mean energy 10 GeV 
in the midst of the LAT energy range. Even in a 10 /iG mag- 
netic field such electrons synchrotron cool quite slowly over 
a timescale of ~ 45 kyr, roughly double the characteristic 
age of the pulsar. Inspection of Figure [3] indicates a spectral 
break at ~ 200 GeV, almost certainly due to a cooling break 
in the electron spectrum. The hard LAT spectrum is therefore 
clearly in the uncooled regime, and so a one-zone model can 
help illuminate this new data. 

We apply a one-zone time dependent spectral en e rgy dis - 
tribution (SED) model, as described in [Abdo et al.l (|2010c). 
This model computes SEDs from evolving electron popula- 



tions over the lifetime of the pulsar in a series of time steps, 
with the energy content of the injected particle population 
varying with time following the pulsar spin down. During 
the free-expansion phase of the PWN (assumed to be ~ 10 4 
years) we adopt an expansion of R oc t, following which the 
radius R oc t 3 , appropriate for a PWN expanding in pressure 
equilibrium with a Sedov phase SNR. Over the pulsar life- 
time the magnetic field B oc 5 , following ~ 500 years of 
constancy. At each time step synchrotron, inverse-Compton 
(Klein-Nishina effects included), and adiabatic losses are cal- 
culated. Synchrotron and IC fluxes are calculated from the 
final electron spectrum. We allow the braking index n of the 
pulsar to vary, thereby changing the age and spin-down be- 
havior of the pulsar. 

We assume the existence of three primary photon fields 
(CMBR, far IR (dust), and starlight) and use th e interstellar 
radiat ion mapcube within the GALPROP suite ( Port er et~aT1 
120051) to estimate the photon fields at the Galactic radius of 
PSR J1826-1334. A distance of 3.9 kpc in the direction of 
the pulsar corresponds to a Galactic radius of 4.7 kpc. At this 
radius, the peak of the SED of dust IR photons corresponds 
to a black body temperature of T ~ 32 K with a density of 
~ 0.9 eVcm -3 , while the SED of stellar photons peaks at 
T ~ 2500 K with a density of - 3.6 eV cm' 3 . 

Spectral measurements consist of LAT and H.E.S.S. data 
points, as well as an estimate of the X-ray spectrum. We 
adopt an X-ray photon index of 2.2 ± 0.3 and a flux of 
(7±2) x 10 -12 erg cm -2 s _1 . The selected ind ex is consistent 
with the indices me asured by XMM-Newt on dGaensler et al.1 
120031) and Suzaku dUchivama et al.l 120091) for the extended 
nebula. The flux lev el is equivalent to sum of all the Suzaku 
regions analyzed by lUchivama et al.l (|2009), with 30 % sys- 
tematic errors assumed. 

A simple exponentially cutoff power-law injection of elec- 
trons, evolved properly over the pulsar lifetime, often pro- 
vides an adequate match to PWNe SEDS. Initially, we fit this 
injection spectrum with four variables: final magnetic field, 
electron high energy cutoff, electron power-law index, and the 
pulsar braking index n. Given the large covariance between 
the braking index and the initial spin period in determining the 
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FIG. 3. — Spectral energy dist ributi on of HESS J1825 — 137 in gamma-rays. The LAT spectral points (in red) are obtained using the maximum likelihood 
method gtlike described in section l4!2l in 6 logarithmically -spac ed energy bins. The statistical errors are shown in red, while the black lines take into account both 
the statistical and systematic errors as discussed in section l4!2l The red solid line presents the result obtained by fitting a power-law to the data in the 1-100 GeV 
energy range using a maximum likelihood fit. A 95 % C.L. upper limit is computed when the statistical significance is lower than 3 a. The H.E.S.S. results are 
represented in blue ( Aharonian et al. 2006). 



age of the pulsar, we fix the initial spin period at 10 ms and 
braking index at 2.5, yielding an age of 26 kyr for the sys- 
tem. This simple injection spectrum slightly underestimates 
the LAT data but the overall fit is still reasonable. For the 
source age of 26 kyr, we require a power-law index of 1.9, 
a cutoff at 57 TeV and a magnetic field of 4 /xG. The corre- 
sponding result is presented in Figure|4](Top). 

Another option to fit the multi-wavelength data is adopting 
the relativistic M axwellian plus pow er-law tail electron spec- 
trum proposed by Spitkovsky (2008). For this injection spec- 
trum, we assume a bulk gamma-factor (70) for the PWN wind 
upstream of the termination shock. At the termination shock 
the ambient pressure balances the wind pressure, fully ther- 
malizing the wind; in this case the downstream post-shock 
flow has 7 = (70 — l)/2. One could also interpret this as 
an effective t emperature kT of m e c 2 (70 — l)/2. Per the 
simulations of Spitkovsky] d2008l) . a power-law tail begins at 
7kT ~ 7/2m e c 2 70, and suffers an exponential cutoff at some 
higher energy. For our modeling we fix the power-law begin- 
ning at ~ 7kT, and allow kT, the power-law index, and the 
exponential cutoff to vary. The best fit, presented in Figure |4] 
(Bottom), is obtained with kT = 0.14 TeV, corresponding to 
an upstream gamma-factor of 5.5 x 10 5 , a magnetic field of 
3 /zG, a cutoff at 150 TeV and a power-law index o f 2.3 close 
to the value of ~ 2.5 proposed bv lSpifkovskvl(l2008l) . The rel- 
ativistic Maxwellian plus power law model matches the multi- 
wavelength data and also directly probes the upstream pulsar 
wind via fitting of 70. 

HESS J 1825— 137 is detected at high significance by the 
Fermi LAT, and demonstrates both morphological similarity 
and flux continuity with the H.E.S.S. regime. The LAT spec- 
tral index of 1.38 ± 0.12 ± 0.16 is consistent with both a 
simple power-law electron injection spectrum, as well as a 
Maxwellian plus power-law injection spectrum over a simple 
power-law. A mean magnetic field of ~ 3 — 4 /iG adequately 
fits the X-ray flux, and an age of ~ 26 kyr is consistent with 
the data. A total of 5 x 10 49 erg injected in the form of elec- 
trons by the pulsar is required to match the gamma-ray flux in 
the nebula. 
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APPENDIX : DESCRIPTION OF SOURCELIKE 

Sourcelike is a tool developed for performing morpholog- 
ical studies of spatially extended Fermi sources. Sourcelik e 
is an extension to pointfit (described in lAbdo et al.l l2010b). 
which was developed to efficiently create Test Statistic maps 
and localize catalog sources with little sacrifice of precision. 
pointfit bins the sky in position and energy and increases ef- 
ficiency by scaling the spatial bin size with energy. Further- 
more, it uses a region of the sky centered on the source whose 
radius is energy dependent: from 15° at 100 MeV to 3.5° at 
50 GeV. It wa s successfully used by and is described in the 
1FGL catalog (Abdo et all 1201 Obi) . Sourcelike deviates from 
pointfit by fitting not the PSF but the PSF convolved with an 
assumed spatial shape. By independently fitting the flux in 
each energy bin, Sourcelike performs an extension analysis 
without biasing the fit by assuming a spectral model. Since 
the PSF ranges a full two orders of magnitude in size over the 
energy range of the instrument, this maximum likelihood ap- 
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FIG. 4. — Spectral energy distribution of HESS Jl 825 — 137 with a simple exponentially cutoff power-law electron spectrum (Top), and a relativistic Maxwellian 
plus power-law electron spectrum (Bottom). The LAT spectral points (red, thin lines denote systematic errors) H.E.S.S. points (blue), and X-ray bowtie (green) 
are shown. The black line denotes the total synchrotron and inverse Compton emission from the nebula. Thin curves indicate the Compton components from 
scattering on the CMB (long-dashed), IR (medium-dashed), and stellar (dotted) photons. 



proach naturally handles the energy dependent PSF and max- 
imizes our sensitivity to extension. 

To find the shape of the source that most closely matches 
the observed photons, the overall likelihood is maximized by 
simultaneously fitting the spatial parameters of the source. At 
each step in the fit, the new shape must be again convolved 
with the PSF. It is for this reason that this optimized software 
was developed. Errors on the fit position and extension are es- 
timated by numerically calculating the curvature of the likeli- 
hood function at the best fit spatial model. The statistical sig- 
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